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Abstract 



We study the influence of Gribov copies, in the Landau gauge, on 
the lattice renormalisation constant Za of the axial current, obtained 
from a Ward identity on quark state correlation functions, with the 
Clover action, in quenched SU(3) gauge theory. A comparison be- 
tween the gauge invariant determination of Za and the gauge depen- 
dent one is discussed. On a 16 3 x 32 lattice at j3 = 6.0 and with 
K = 0.1425, the values, on a sample of 36 configurations, are: Za 
= 1.08(5) (gauge dependent calculation) and Za = 1.06(2) (gauge 
independent calculation). We find that the residual gauge freedom 
associated to Gribov copies induces observable effects, which, at the 
level of numerical accuracy of our simulation, are included in the sta- 
tistical uncertainty inherent in a Monte Carlo simulation. Doubling 
the statistics suggests that the fluctuation due to the lattice Gribov 
ambiguity scales down at least as fast as a pure statistical error. 
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1 Introduction 



In continuum non Abelian field theories, most popular choices of fixing the 
gauge (e.g. Landau, Coulomb) suffer from the Gribov ambiguity fl]]. It is 
now well established that this problem also affects the lattice formulation of 
these theories [[J-0 . This problem has been neglected for a long time be- 
cause, in principle, the computation of gauge invariant operators in compact 
lattice theories does not require gauge fixing. Fixing the gauge is, however, 
necessary in several cases. Monopole studies in SU(2) pure gauge theory 
have been done in the unitary gauge and the effect of the Gribov ambiguity 
on the number of SU(2) monopoles has been investigated jS], [(J. The authors 
conclude that, in their case, the Gribov noise does not exceed the statisti- 
cal uncertainty. In SU(3) gauge theory, gauge fixing is essential in the the 
computation of gauge dependent quantities, such as gluon and quark prop- 
agators. There are now several studies of lattice propagators. The gluon 
propagator has been calculated in 0-0 with the aim of studying the mech- 
anism through which the gluon may become massive at long distances. More 
recent attempts have investigated its behaviour as a function of momentum 



IffC] , [TTH . Analogous studies have also been performed on the quark propa- 
gator (see, for example ||). In practice, there are also cases in which it is 
convenient to implement a gauge dependent procedure for the computation 
of gauge invariant quantities [12]j-||14||. For example, smeared fermionic in- 



terpolating operators are widely being used in lattice QCD spectroscopy and 
phenomenology, in order to optimise the overlap of the lower-lying physical 
state with the operator. The point-splitted smeared operators are gauge de- 
pendent, and therefore the gauge must be fixed before they are calculated. 
In particular, the calculation of the decay constant of the B meson in the 
static approximation, in which the 6-quark has infinite mass, requires the 
computation of the two point correlation function of the axial current. The 
isolation of the lightest state at large times is not possible if local (gauge 
invariant) operators are used. A nice way out consists in smearing the bilo- 
cal operator over a small cube and extracting fg by forming suitable ratios 
of smeared and local correlation functions ||13|| . This is an explicitly gauge 
dependent procedure which is most naturally carried out in the Coulomb 
gauge. In ref.JT^] the smeared - smeared correlation functions on a few indi- 
vidual configurations were computed. Two Gribov copies were produced per 
configuration. The Gribov noise on individual configurations was found to 
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vary from 1% to 50% depending on the time-slice, which implies that it may 
still be a considerable effect after averaging over configurations. However, it 
was not possible to estimate its effect beyond individual configurations. The 
reason is that in such a study other sources of error dominate, such as the 
systematic error arising from fitting the exponential decay of the correlation 
function with time. Thus the isolation of the Gribov noise is difficult. [] 

In this paper we study a different physical quantity, namely the renormal- 
isation constant Za of the lattice axial current. A knowledge of these renor- 
malisation constants is necessary for matching the matrix elements computed 
using lattice simulations to those required in a definite continuum renormal- 
isation scheme. Provided that the lattice spacing is sufficiently small it is 
possible to calculate these renormalisation constants in perturbation theory. 
For a more reliable determination of these constants it has been suggested to 
impose the chiral Ward identities of QCD non-perturbatively O, |TS|. Here 
we focus our attention on the determination of the role of the Gribov ambigu- 
ity in the calculation of Za, obtained from quark state correlation functions. 
A recently proposed method to determine Za and other renormalisation con- 



stants, based on truncated quark Green functions in momentum space |fL9 
can also in principle be afflicted by Gribov fluctuations. 

Since reasonably small errors are expected, in this kind of calculations, it 
is crucial to investigate the role of the Gribov noise. Moreover, the renormali- 
sation constant Za of the axial current is particularly well suited to the study 
of the Gribov fluctuations, mainly for two reasons. Firstly, Za can be ob- 
tained from chiral Ward identities in two distinct ways: a gauge independent 
one, which consists in taking the matrix elements between hadronic states, 
and a gauge dependent one, which consists in taking the matrix elements 
between quark states. Hence, there is an explicitly gauge invariant estimate 
of Za which is free of Gribov noise and which can be directly compared to 
the gauge dependent, Gribov affected, estimate. The second advantage is 
that Za is obtained by solving a first degree algebraic equation for each lat- 
tice time slice, thus avoiding the usual systematic errors arising from fitting 
exponentially decaying signals in time. 



1 More recently, a high statistics study of Jb in the static limit |16 uses a different 
method for constructing ratios of smeared and local correlators which avoids fitting. This 
method, however, requires a large temporal extention of the lattice. 
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2 Chiral Ward identities for Za 



The theoretical framework for the non-perturbative evaluation of Za for Wil- 
son fermions, has been developed in [18|]. The renormalisation constant is 
obtained through Ward identities generated by axial transformations. A 
first application of these techniques in numerical simulations using the Wil- 
son action was attempted in P0fl . The extension of these methods to the 
0(a) improved Clover action (a is the lattice spacing) was presented in 
which we follow most closely. Here we only give a brief outline of the 
results which are essential to our work. 

In this study terms that, close to the continuum limit, are effectively of 
0(a) are eliminated by using the Clover action [^3j and rotating all quark 
fields of the matrix elements according to the "improved improvement" pre- 



scription of [24 



va — * — — va < — 

i> — (1 - y7 ■ D)tp ; V — i>(l + y7 ■ D) 



(1) 



(r is the Wilson term parameter; in this work r = 1). The two fermion local 
operators considered in the following are the axial and vector currents and 
the pseudoscalar density: 

f - * f 

f - ^ 



P((x) = i)(x)^'—i)(x) 



(2) 



(/ is a flavour label and the notation is generic for any quark fields ip and 
ip). In order to ensure that the lattice axial current has the correct chiral 



properties, it is normalised by a renormalisation constant Za ||18|| ; this implies 
that 



2p 



d 4 fd 3 y< A{(y,t y 



)P£ f (0,0) 



> 



fd*y< Pi(y,t y )P^(6,0) > 



if i 



2m 

~Z~A 



(3) 



(m is the bare quark mass and \P > the pseudoscalar state). As can be seen 
from the above equation, p is gauge invariant. 
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The gauge dependent calculation of Za relies on the following Ward iden- 
tity for quark Green functions [18, 2^ 



2pTr 



— pra ) Tr 



d 4 x / d 3 y < u a (y) (u(x)^d(x)) dp(0) > 



d 3 y < (j 5 d(y)d(0) + u(y)u(0)j 5 ) > afi 



(4) 



In eq.(|4]) we work explicitly with up and down quark fields with spinor labels 
a and (3. The trace is over colour indices. The expectation values on both 
sides of (|J) are evaluated as functions of t y . Taking the value of 2p obtained 
using eq.(|3]), Za can then be determined. In order to enhance the signal, 
we add in both sides of (|]) the four contributions (a, /3)=(1,3), (3,1), (2,4) 
and (4,2), which were found to give the clearest signal p2| . A plateau in t y 
is typically obtained and Za is estimated from it. The crucial point is that 
both sides of eq.(||) are gauge dependent, and thus this determination of Za 
is in principle sensitive to the Gribov noise. 

A gauge invariant determination of Za is obtained through the Ward 
identity 

2p J d A x J d 3 y< Pi(x)Al(y)V p h (0) > 

A -pra\f f91 Jd 3 y< Vl(y)V p h (0) > 
1 




-'I — - pra) f hl I d 3 y< Al(y)A l p (Q) > 



(5) 



In the above equation, the vector current renormalisation constant Zy is also 
needed. For the Clover action, Zy is calculated with the aid of the so-called 
conserved and improved vector current 



1 rv 



^(x)(7 M - r)t^(x)^(x + p) + ip(x + /})(7„ + r)Ul(x)tjj(x) 



+ {x 



f ✓ _ 



(6) 



Since the current is conserved, its renormalisation constant is precisely 1, and, 
since it is "improved", its matrix elements have no corrections of 0(a). The 
normalisation constant Zy of the local vector current VJx) is determined 
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through the ratio of the vacuum to vector state matrix elements of the two 
vector currents f25| : 

„ _<oK OI (o)\v»> 



' v 



(here \ V^ > denotes the vector state). By calculating all correlation functions 
of eq.@ and Zy from eq.(|7]) and by requiring that eq.® holds at each t y , 
we can determine Za in an explicitly gauge invariant fashion for which the 
Gribov ambiguity is irrelevant. 

One more comment is in place here: the terms proportional to pra on 
the right- hand- side of eqs. (f|) and (|5|) arise from the rotations of the fermion 
fields defined in eq. ([!]) , which are inherent to Clover action improvement [p3 



The 7 • D rotations, combined with the equations of motion, generate contact 



terms which, to 0(a), give rise to the terms proportional to pra [p2 



3 Lattice Gauge Fixing and Gribov Copies 

Gauge fixing and the generation of Gribov copies on the lattice is by now a 
standard procedure. Given a thermalised configuration generated by a Monte 



Carlo simulation, the Landau gauge is fixed by minimising the functional [26 



F[17»] = — ^JfeTrSEW ( g ) 

with respect to g. Here V is the lattice volume and U 9 (n) = g(n)U fJi (n)g(n + 
pY is the compact SU(3) gauge field, gauge transformed by a local gauge 
transformation g(n). The extrema of F correspond to configurations that 
satisfy the gauge condition d^A 9 ^ = in discretised form. The minimisa- 
tion of F is obtained through iteration: each lattice site is visited and F is 
minimized locally. After several lattice sweeps F becomes constant, and the 
gauge is fixed. This is nothing else but a discretised analogue of the contin- 
uum formulation [JS7J, according to which the local minima of the functional 
jFfyl^] = — Ylt=i Tr J d^xi^A 9 ^) 2 with respect to g(x) correspond to configura- 
tions in the Landau gauge. Even if there is a straightforward correspondence 
between the lattice and continuum gauge fixing procedures, it is interest- 
ing to note that on the lattice, because of discretisation, there can be more 



minima than in the continuum [28 
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The production of Gribov copies consists in generating gauge equivalent 
configurations, by applying random gauge transformations to our thermalised 
link configuration |2B|, [Q. The gauge is then fixed and F is measured. 
Since it is a gauge dependent quantity, a different value of F for two gauge 
equivalent, gauge fixed configurations means that they are Gribov copies. 
The over-relaxation technique of (which consists in accelerating the gauge 
fixing algorithm by raising the gauge transformation g(x) to a real tunable 
power u at every iteration) has been implemented at fixed uo. The over- 
relaxation itself can be used for the generation of Gribov copies, by varying 
the value of u, as proposed in In this work we have opted for the random 
gauge transformation method. 



4 Results 

We work in the framework of the quenched approximation with the Clover 
action of SU(3) gauge theory. The lattice volume is V = 16 3 x 32 and 
(3 = 6.0. After 3000 thermalising sweeps, 36 configurations were generated, 
separated by 1000 sweeps. An 8 hit Metropolis algorithm was used. For 
each thermalised configuration, we generated 6 Gribov copies. This was 
done by fixing the gauge both on the original configuration and on 5 gauge 
equivalent replicas, obtained by applying random gauge transformations. It is 
remarkable that, in our case, each random gauge transformation produced a 
Gribov copy. This high probability to find Gribov copies is a characteristic of 



large volume lattices in the confined region, as discussed in [|30| , [3T|| . We fix the 
Landau gauge, using the over-relaxation algorithm suggested in [F29] at fourth 
order in the over- relaxation parameter to = 1.72. The stopping condition we 
have imposed is that the relative variation SF/F, of the minimised functional 
F, between two consecutive gauge fixing sweeps be less than 10 -8 . This value 
guarantees a good quality of the gauge fixing, allows us to distinguish Gribov 
copies, and it is typically reached after a number of sweeps which varies 
between 500 and 1500. The gauge fixing was done on an IBM Rise 6000/550 
equipped with 128 Mbyte of RAM memory and with a CPU working at 
42 MHz; with this machine a single gauge fixing sweep takes about 40 s. 
The quark propagators, rotated as indicated by eq.fll]) were obtained at a 
Wilson hopping parameter value of K — 0.1425, which, for the Clover action, 
corresponds to a pion of roughly 900 MeV. 
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Before passing to a detailed discussion of the Gribov noise, it is appro- 
priate to present a comparison of the results for Za, obtained in a gauge 
invariant way (see eq.(|5])), to those based on the gauge dependent identity 
of eq.(||). This calculation has been already presented in p2f , on the first 18 
configurations of our ensemble; here we have doubled the statistics. More- 
over, we have re-gauged the first 18 configurations in order to reach the 
better quality of gauge fixing required for this study. In Fig.(l) we show 
the behaviour of the two estimates of function of t y , calculated on 

the same ensemble of 36 configurations. The gauge invariant values of Za 
show a flat behaviour already at t y = 5 and with small error bars. The new 
value of Za, obtained with the gauge independent technique, over 36 config- 
urations, is Za = 1.06(2) to be compared with the old one obtained over 18 
configurations: Za = 1.09(3) (see ref . p2| ) . 



In the gauge dependent case, we have taken the average over the Gribov 
copies, in the way that will be discussed below. In this case, the Za behaviour 
becomes flat only at t y = 10 showing a large sensitivity to the contact terms 
of the Ward identity at small t y values. The new value of Za, obtained 
from the gauge dependent Ward identity, as can be seen from Fig.(l) is 
Za = 1.08(5); to be compared to the value obtained from 18 configurations, 
Za = 1.14(8) (see ref.f22|). As already stressed in f22|, this method gives 
results that are compatible within the errors with the gauge independent 
ones. Moreover, the error of the gauge invariant calculation of Za is always 
smaller than the error of its gauge dependent counterpart. This is due to 
the fact that the quark state correlation functions fluctuate more than the 
gauge invariant correlation functions, but it may also indicate the presence of 
another effect, which is probably the Gribov ambiguity. In the hypothetical 
case of two determinations of Za, affected by the same statistical error, the 
gauge dependent estimate should fluctuate more due to the Gribov noise. 
Then the amount of Gribov noise could be estimated as the difference (in 
quadrature) of the two errors. 

Normally, in a standard simulation of gauge dependent quantities, one 
does not generate Gribov copies. Consequently, one measures a given quan- 
tity by taking the average and error over the ensemble of the gauge fixed 
configurations that have been generated. This implies a particular and arbi- 
trary choice of Gribov copies. The error estimated, for example, by a jacknife 
method, is not purely statistical as it implicitly contains the uncertainty due 
to the particular choice of a Gribov copy. 
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In our case, having generated G = 6 copies for each of the N = 36 
thermalised configurations, there are G N possible combinations that we may 
consider when forming a particular ensemble. To the best of our knowledge, 
the distribution of the Gribov copies of a given configuration is unknown; 
thus the weight to be associated to it is arbitrary. Moreover, any technique 
used to generate different Gribov copies selects a particular copy in a com- 
pletely uncontrolled way. Hence, we assume that the particular choice of 
different combinations of Gribov copies when forming a statistical ensemble 
is arbitrary. In order to exhibit the effect of such arbitrariness, we show in 
Fig. (2) the behaviour of ZjS$>y) for 4 arbitrary choices of copies. The 4 differ- 
ent behaviours are compatible, and the same is true for the jacknife errors. 
It is clear, however, that the presence of Gribov copies is a visible effect; each 
of the six estimates of Z& shown has a slightly different profile as a function 

of ty. 

We now expose the procedure we implemented for taking into account the 
Gribov ambiguity. The gauge dependent traces of the two and three point 
correlation functions appearing in eq.(|]), calculated on a single configuration 
c and for a particular Gribov copy g are denoted by: 



(in the above equations the Dirac indices have been implicitly averaged over, 
as explained in Sect. (2)). Then, for a given configuration c, we consider our 
"best estimate" of these matrix elements to be their average over the G = 6 
Gribov copies: 



The average of the above T's over the 36 configurations will be taken as our 
"best estimate" < T > for the gauge dependent traces. On the other hand, 
p, being gauge invariant, does not depend on a particular choice of Gribov 
copies, but only on the configuration ensemble. Then Z^(t y ) is obtained by 
applying eq.(|3]) as follows 




(9) 



1 



G 




(10) 



Z A{ty) 



2p(t y ) 



< T 3 (ty) > 

< T 2 (ty) > 



+ p{t y )ra 



(11) 
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The error is obtained by a standard jacknife method performed on the quan- 
tities < T 2 (t y ) >,< T 3 (t y ) > and pity), by decimating one configuration at a 
time. This is the gauge dependent ^(iy) result shown in Fig.(l). 

We want to stress that the values of < T 2 > and < T 3 > fluctuate 
more than their ratio < T 3 > / < T 2 >. The latter quantity is gauge 
invariant, according to eq. (|TI|) . For example, at t y = 10, 5T 2 / < T 2 >= 5.6%, 
8f 3 / < T 3 >= 5.1% and < g> 3 /<^ > = 2.3%. The strong reduction of the 
relative error indicates a great sensitivity of T 2)3 to the Gribov fluctuation, 
as opposed to a relative stability of the gauge invariant quantities. 

In order to estimate the uncertainty arising from a particular choice of 
copies, out of the 6 36 possible ones, we have applied a procedure which takes 
this arbitrariness into account. We have chosen randomly 10 4 combinations 
of copies of our ensemble and have calculated Z^ity) for each one of these 
combinations at fixed t y . The histogram of the values obtained for ZA{t y ) is 
well fitted by a Gaussian, the r.m.s. width of which is taken as an estimate 
of the fluctuation. In Fig. (3 a) we compare the jacknife error of our "best 
estimate" to the width of the Gaussian. We see that for all t y the width is 
smaller than the jacknife error. This implies that the fluctuations induced 
by the particular choice of Gribov copy when forming the ensemble are small 
and do not overcome the statistical uncertainty. 

As it is also important to understand how the above behaviour scales 
with increasing the number of configurations, we have performed the same 
analysis for the first 18 configurations (half of our ensemble). The result 
is shown in Fig. (3 b). Comparing Fig. (3 a) to Fig. (3 b), we note that both 
errors scale at least as \^2. Thus, within our moderate statistics, we find that 
the error, even if affected by the Gribov ambiguity, decreases with increasing 
configuration number. 

In conclusion, our investigation, even within its limitations, shows that, 
on the lattice, the residual gauge dependence associated with the Gribov 
copies does not generate an overwhelming fluactuation of the Za measure- 
ments performed by the gauge dependent method. Nevertheless the arbi- 
trariness associated to a particular choice of Gribov copies is a visible effect 
which manifests itself, especially in the behaviour of function of £„. 

The jacknife error is greater in the gauge dependent determination than in 
the gauge independent one. However, in the former case, even though the er- 
ror is afflicted by the Gribov uncertainty, it is still decreasing with increasing 
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statistics. This implies that the uncertainty arising from the Gribov am- 
biguity may be a secondary effect. Analogous studies on different physical 
quantities and renormalisation constants could further support this belief. 
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6 Figure Captions 

FIGURE 1: Comparison of the gauge independent calculation of Za^v) (di- 
amonds) to the gauge dependent one (crosses). The errors are jacknife. The 
crosses have been slightly displaced to help the eye. 

FIGURE 2: The gauge dependent calculation of ^a(^) for 4 arbitrarily cho- 
sen Gribov copies. The errors are jacknife. 

FIGURE 3: Comparison of the jacknife error (crosses) to the r.m.s. Gaussian 
width (diamonds) due to the Gribov ambiguity (see text). The crosses have 
been slightly displaced to help the eye. Case (a) is with 36 configurations; 
case (b) with the first 18 configurations 
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